Towards HPC simulations of billion-cell reservoirs by multiscale mixed methods

A three dimensional parallel implementation of Multiscale Mixed Methods based on non-overlapping domain decomposition techniques is proposed for multi-core computers and its computational performance assessed by means numerical experiments. As a prototypical method, from which many others can be derived, the Robin Coupled Method chosen explained in detail. Numerical results problems ranging millions up to more than 2 billion cells highly heterogeneous anisotropic rock formations SPE10 benchmark are shown. The relies direct solvers both local interface coupling system. We find good weak strong scalalability as compared against state-of-the-art global fine grid solver Algebraic Multigrid preconditioning single two-phase flow problems.